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We investigate a model of a two-component Bose-Einstein condensate residing on an optical lattice. Within a 
Bogolioubov-approach at the mean-field level, we derive exact analytical expressions for the excitation spectrum 
of the two-component condensate when taking into account hopping and interactions between arbitrary sites. 
Our results thus constitute a basis for works that seek to clarify the effects of higher-order interactions in the 
system. We investigate the excitation spectrum and the two branches of superfluid velocity in more detail for 
two limiting cases of particular relevance. Moreover, we relate the hopping and interaction parameters in the 
effective Bose-Hubbard model to microscopic parameters in the system, such as the laserlight wavelength and 
atomic masses of the components in the condensate. These results are then used to calculate analytically and 
numerically the drag coefficient between the components of the condensate. We find that the drag is most 
effective close to the symmetric case of equal masses between the components, regardless of the strength of the 
intercomponent interaction and the lattice well depth. 
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I. INTRODUCTION 



The emergence of laser cooling techniques and their appli- 
cations to realizing the phenomenon of Bose-Einstein conden- 
sation (BEC) in the laboratory, has paved the way for a study 
of the rich physics present when atoms condense at ultralow 
temperatures on an optical lattice [T][2][3]. The BEC itself is 
a coherent matter wave, and has attracted much attention both 
theoretically and experimentally over the past decade. One of 
the remarkable features of a BEC residing on an optical lattice 
is the extent to which physical quantities such as tunnel cou- 
pling and on-site interaction may be controlled experimentally 
simply by adjusting the lattice parameters. This is done by 
controlling the interference pattern of the lasers setting up the 
optical lattice. For instance, by causing the depth of the lattice 
potential to increase, when would expect a resultant decrease 
of the hopping amplitudes and an increase of the on-site inter- 
action. The possibility to alter the lattice parameters directly 
during the experiment, and thus influencing the physics, is 
clearly intriguing. Moreover, experiments carried out on such 
systems are extremely well controlled since there is no disor- 
der present. Since the atoms reside on an enginereed lattice, it 
is possible to investigate the physics by means of standard the- 
ories in condensed matter physics, such as the Bose-Hubbard 
model [4j. As pointed out in Ref. [3|, BECs residing on op- 
tical lattices have several advantages compared to ultracold 
atoms in a non-condensed phase. The main point is that the 
temperatures and densities for ultracold atoms and BECs both 
differ by three to four orders of magnitude. One consequence 
of the much higher particle densities for BECs is that atomic 
interactions become crucial with regard to the physics. 

By allowing for more than one component of bosonic atoms 
on an optical lattice, one opens up an exciting avenue of 
physics to explore J5] |6] [7] [8] |9] [10) . The physical realization 
of such a multicomponent BEC includes condensates with 
spin degrees of freedom (spinor condensates), two or more 
hyperfine states of the same atomic species that condense si- 



multaneously, or simply two distinct atomic species. The two- 
component condensate has been shown to be a more rich en- 
vironment to explore than a single-component BEC due to the 
possibility of an "entrainment" coupling between the conden- 
sate components, see Ref. |9l[T0]and references therein. Such 
a system may be studied at a mean-field level by employing a 
Bogolioubov-approach, which may provide information about 
both the transition from a superfluid to Mott insulating state 
and also the quasiparticle excitation energies which arise from 
the condensate. By means of the Landau criterion, it is also 
possible to obtain information about the superfluid velocity of 
the condensate from the excitation spectrum. 

Very recently, the excitation spectrum for a two-component 
Bose-Einstein condensate was obtained for a limiting case 
in Ref. (12). In that work, the author presented a correc- 
tion to erroneous results previously reported in the literature 
lfT3ll . The calculations were performed under the standard 
assumptions of nearest-neighbor hopping and on-site interac- 
tions only. It would clearly be of interest extend calculations 
beyond these approximations, in order to investigate how the 
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FIG. 1: (Color online) An optical lattice setup by counter- 
propagating lasers serves as a potential landscape for two atomic 
species, denoted by the red (dark) and green (light) spheres. Each 
species of atoms may hop from site to site and also interact with both 
inter- and intraspecies atoms. 
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excitation spectrum is affected by taking into account longer- 
range hopping and longer-range interactions. One of the pur- 
poses of the present paper is to extend the calculations of Ref. 
\V2\in this direction. 

Another goal in this paper is to address the effect of drag 
between the atomic components in a two-component Bose- 
Einstein condensate residing on an optical lattice. Such a drag 
effect points to a mutual transfer of motion between the com- 
ponents, and was first investigated in 3 He- 4 He superfiuid mix- 
tures by Andreev and Bashkin |1 1]. In Ref. |7|. the drag effect 
for a two-component Bose gas was explored in the continuum 
limit. We will here derive an analytical expression for the in- 
tercomponent drag in a two-component Bose-Einstein con- 
densate residing on an optical lattice, and relating it directly to 
the microscopic parameters in the system which are possible 
to tune experimentally. 

We organize this paper as follows. In Sec. [TTJ we establish 
the theoretical framework to be used in deriving our main re- 
sults. In Sec. Ill we provide an analytical solution for the 



excitation spectrum of a two-component Bose-Einstein con- 
densate for arbitrary hopping and interaction between sites 



(Sec. Ill A i and investigate the superfiuid velocity and phase- 



separation condition in more detail for two limiting cases in 



Sec. Ill B and III C Also, we present a correction to the con- 
dition for phase-stability of the two components, which de- 
termines whether the species are spatially miscible or not. In 
Sec. 



HID 



we first relate analytically the parameters in the 
two-component Bose-Hubbard model directly to the funda- 
mental physical quantities such as mass and trapping poten- 
tial. Then, we combine these results with the expressions for 



the excitation energies in Sec. IIIB and obtain an analytical 



equation for the drag coefficient in the system. The drag coef- 
ficient is then studied as a function of the microscopic param- 
eters. Finally, we give concluding remarks in Sec. |IV] The 
system under consideration is shown schematically in Fig. [T] 



H. THEORY 

The starting point for our calculations is a microscopic 
Hamiltonian H for an ensemble of bosonic atoms that are con- 
fined by a slowly varying external harmonic trapping potential 
Vj t a(r) and subject to an additional optical lattice potential 
Vo(r). In terms of boson field operators ip a (r), where a de- 
notes the boson-component, H may be written as (K = 1) 



2m a 



H = J2 [ drVl(r)[- 
x V-aO) + - [ dr "tPl( r )^l( r h a i3iPi3(r)ip a (r), (1) 



a/3 



Here, a a , clab are intraspecies and interspecies s-wave scat- 
tering lengths. The interaction strength is assumed to be re- 
pulsive and, in general, different for each of the boson compo- 
nents: {a a ,aAB} > 0. To obtain a second-quantized Hamil- 
tonian in a lattice-formulation, we assume that the field op- 
erators vp a (r) may be expanded in a Wannier function ba- 
sis set. The physical motivation for this is that the bosons 
are assumed to spend most of their time in the minima of 
the optical lattice potential, with occasional tunneling from 
one site to another. In this case, a set of localized Wan- 
nier functions where only the lowest lying excitation level is 
taken into account is expected to be a reasonable choice of 
basis. We consider here a two-dimensional model, such that 
''Pair) = J2i b iaW a (x - Xi)w a (y - yi). Here, b ia are bo- 
son annihilation operators for species a on the lattice point 
i, while W a (r — r,) = w a (x — Xi)w a (y — yi) are single- 
particle Wannier states for boson species a centred around 
lattice point i at (xi,yi). Inserting this expansion into Eq. 
([T]l yields an effective Bose-Hubbard like model, defined by 
the Hamiltonian 



a i^j ia 

+ 2 ^ 5^ Uijklapblj^pbkpha- 
ijkl a (3 

The parameters of this model are expressed as 

V 2 



(2) 



t ija = ~ / drW*ir-n) 



+ V (r) 



drW*(r-r<) 



2m a 

+ V T , a (r)^W a (r-r 3 ) 1 

- V 2 
- o ^« + Voir) 

+ Vt 1 „(r)lw a (r-r i ), 



Uijkiap = l at 3 J drW*ir - r,)Wpir - r,) 

x W P ir-r k )W a ir-n). (3) 

So far, we have made no approximations apart from the 
assumed field expansion. The integrals given above may 
be evaluated analytically by specifying the explicit form of 
W a ir). Let us consider the following generic form for the 
trap and laser potential: 



y T , Q (r) = ^(^ 2 +c,y + ^z 2 ), 
V (r) = V x sin 2 (k x x) + V y sin 2 (k y y) + V z sin 2 (k z z). 



(4) 



where jap denotes the onsite-interaction for both boson 
species, m a is the mass of boson species a, and fi a is its 
chemical potential. Specifically, we have Q 



la/3 



Aira a /m a ifa = /3 

2n(mA + m B )aAB/niAiTiB if a ^ (3 



Here, ojj is the frequency of harmonic trapping potential as- 
sociated with the j-direction while the wave vector kj for the 
optical lattice is related to the wavelength A of the laser light as 
kj = 2ir/\j, such that the lattice period becomes aj = Aj/2, 
j G {x,y,z}. In the harmonic approximation (4j[T7], where 
the bosons have a small probability of being located far from 
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each lattice site and higher energy states in each lattice poten- 
tial may be neglected, the exact Wannier functions can be re- 
placed with their harmonic-oscillator approximation to a sat- 
isfactory degree. Then, one may write 



v a (x - Xi) = y - J 



1/4 



Wi, a = \Ju 2 x + 2V x k%/m a . 



(5) 



and similarly for w a (y — y.;) and w a (z — Zi). Since the Wan- 
nier functions are known, one may derive analytical expres- 
sions that relate the parameters in Eq. Q to the microscopic 
parameters in the system. For the hopping term tij >a , previous 
works have neglected the influence of the trapping potential on 
this parameter by demanding that Vx,a (v) varies much more 
slowly than Vo(r). In this work, we derive a more general 
expression for both the hopping parameter and the interaction 
term by generalizing previous results to the two-component 
case and also by including the effect of the trapping potential. 
This is done towards the end of Sec. [In] 



III. RESULTS 

We now proceed to derive an analytical expression for the 
excitation energies of the elementary quasiparticles of the con- 
densate. The standard approximation consists of only consid- 
ering nearest-neighbor hopping and on-site interactions. To 
begin with, we include all orders of hopping and interactions 
without any site-limitation. We then explicitly consider two 
cases of particular relevance. Finally, we relate the micro- 
scopic parameters of the system to the hopping and interaction 
term in the effective Bose-Hubbard Hamiltonian. 



A. General solution 

By introducing a mean-field decomposition of the interac- 
tion terms allows us to consider the case where a macroscopic 
number of particles have condensed into the zero-momentum 
state. Let us define the Fourier-transformed boson operators 



h \ , U p—ikri 

u ia prj / u ka K ? 



(6) 



which inserted into Eq. Q may be written as 

H = J> fc ,a + T a )b{ Q b ka + ^Y1 [Z) 2^ fe2 ' ks > ki) 



hex 



where we have defined the generalized intraspecies potential 
U a (k 2 ,k 3 , fc 4 ) = U a (0, 0, 0) + U a (Si,S 2 , S 3 ) 

{Si} 

x gi(fc2-<5i — fc 3 <5 2 — fc<r<53) (g) 



and the interspecies potential 

U A B(k u k 2 , k 3 , fe 4 ) = U AB (0, 0, 0) + U AB (S 1 ,S 2 ,S 3 ) 

{<*»} 

x ^ e i(k 1 -S 1 -k2-S 2 -k4-S 3 ) _|_ e i(fc3-<5i-fc4'<52-fc2'<53)^ (g^ 

Above, the quantities U a (8i, 82,83) and Uab(8i, 82, 83) de- 
note the interaction strengths and their dependence on the site 
distance between the particles involved in the scattering pro- 
cess, while N s denotes the number of lattice sites. Also, we 
have assumed that the energy off-set at each lattice site is sim- 
ply a constant ei a = T a . The interactions are related to the 
scattering potential Uijt-iap as follows 

U a (8i, S 2 , 83) — Ui^i+S 1 ,i+S 2 ,i+S 3 ,aen 
Uab(8i, 8-2,83) — Ui^ + S ly i+S2,i+S 3 ,AB, 

(10) 

and are thus assumed to be independent on at which particular 
lattice site i the scattering takes place, as is reasonable. The 
kinetic energy term is given by 



s 



t a (8)e 



-ikS 



(11) 



where the summation over 8 is to be taken over all neigh- 
bor sites. In Eqs. ^ and ([9J, the summation over {8i} = 
(81,82,83) is to be taken over all possible combinations of 
on-site and off-site lattice points except for pure on-site scat- 
tering { 8i} — {0}. In this way, the first term in Eqs. ([8| and 
(|9]l represents the on-site interaction while the second term 
incorporates scattering involving multiple sites. 

Since we are considering the condensed phase, we may 
write 

&0afyL = 4q^0q + 1 - N 0a > 1, 

N 0a = N a -J2b{ a b ka , (12) 
fe 

where the ' superscript over the sum denotes summation over 
all modes except fe = 0. Physically, we are stating that the 
number of atoms in the zero-mode state k — dominates the 
contribution to the total number of atoms for all fc-modes. The 
biquadratic terms may be reduced to bilinear form by retaining 
only the interaction between the fe = modes and other fe ^ 
modes. Since the number of atoms No a in the fc = mode for 



x KtcP^cfikzabk^ki+^M+ki + U A B(k 1 ,k 2 ,k3,k 4 ) b< 
x bl lA bk 2 Ab k3B b kiB < ; >k 1 +k 3 ,k2+k4 1 ( 7 ) 



atom species a is assumed to satisfy Eq. ( 12 1, we may replace 



Next, we explicitly take into account the <5-function con- 
straints on the particle momenta in Eq. which allows us to 
reduce the Hamiltonian to a sum over the atom species a and 
a single sum over momentum fe. In this way, one obtains 



(13) 
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where we have defined 



and the interaction terms 



r N 2 
H = Y, l N <x(T a + e ,a) + ^-^(0,0,0) 



^^(7^(0,0,0,0) 



(14) 



H a = el + (n a /2)[U a (k, 0, fc) + U a (0, fc, 0) + U a (0, 0, fc) + U a {k, fc, 0) - 20* (0, 0, 0)}b{ a b ka 
+ (n a /2)[U a (-k,0,0)bl a bl ka + U a (0,k, -k)b 



(15) 



Hab = [U AB (k, 0, 0, k)b{ A b kB + U AB (0, fc, fc, 0)b kA b{ B + U AB (fc, 0, -fc, 0)6^^ + 0^(0, fc, -fc, 0)b kA b- kB ] 

x \fnAJiB + n A [U AB (0,0, fc, fc) — 17.45(0, 0, 0, 0)]b kB b kB + n B [Uab (fc, fc, 0, 0) - U AB (0, 0, 0, 0)]^ A 6 feA , (16) 

where e£ = £fc Q + J^s The above equation describes the Hamiltonian of a two-component Bose-Einstein condensate 

residing on an optical lattice with a drag between the atomic species. By diagonalizing Eq. ( 13 1, we obtain the quasiparticle 



spectrum which allows for a further study of the different phases that may be expected for the condensate and also how the 
superfluid velocity depends on the interaction parameters. Using the basis 

4>k = [bkA, b-kA, b kBl b- k B,b hA , bl kA , bl B ,b^ kB } T , 



the Hamiltonian can now be written in compact matrix form: 



(17) 



(18) 



where the matrix M. k reads 



upon defining the auxiliary matrices 



Mi(fc) = 



M 2 (fc) = 



(E A {k) 


V o 
/ 

17a (fc) 


\v 2 *(k) 



Mi(fc) M 2 (fc) 
M 2 (fc)* Mi(fc)* 



Vi(k) 

E A (-k) Vi(-fc) 

E B (k) 

Vt{-k) E B (-k) y 

U A {k) o vj(ky 

V 2 *(~k) 

V 2 *{-k) U B (k) 

U B (k) 



We have introduced the following notation: 



E A (k) 



E B (k) = 4 



U A rr-7 



[U A (k, 0, fc) + U A (0, k, 0) + U A (0, 0, fc) + U A (k, fc, 0) - 2U A (0, 0, 0)] 
n B [U AB (k, fc, 0, 0) - U AB (0, 0, 0, 0)], 

~[U B (k, 0, fc) + U B (0, fc, 0) + U B (0, 0, fc) + U B (k, fc, 0) - 2U B (0, 0, 0)] 



(19) 



(20) 



+ n A [U AB (0, 0, fc, fc) - U AB (0, 0, 0, 0)], 
Uj(k) = n,jUj(—k, 0, 0) with j = A, B, Vi(k) = ^n^nBU AB {k,Q,Q,k), V 2 (k) = ^/n^U AB {Q, fc, -fc,0). (21) 



In order to obtain Eqs. ( 19 1 and (20 1, we made use of the fact that the matrix M. k a^ must be Hermitian, since the eigenvalues 
have to be real (see discussion below). Our ultimate goal is to obtain a Hamiltonian that may be written as 



H=H a + -Y,<S>lVk<S>k, 



(22) 
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where the matrix T>k contains the excitation energies. Note that Hq will in general be different from Ho. The new basis 
is related to the old one through the diagonalization matrix Tk, and also satsifies the correct boson commutation relation: 
® k = T},4>k, ^k&\, — {^k^k) T = <?3- From the requirement that the new basis also consists of boson operators, one finds that 
the relation Ttasfk — 03 must be satisfied. From this, one may infer that (M.k<iz) — fk(DkO'3)T I ^ 1 , which means that Tk 
diagonalizes the matrix (A^feda). The corresponding eigenvalues are contained in the matrix L>k&3, and may be determined by 
considering \Mk&3 — Al| = 0. Evaluating the above determinant yields four distinct eigenvalues = ±£fc CT , a = ±1. 

Before carrying out the diagonalization procedure, it is advantadgeous to make a simplifying observation: if the interaction 
potential satisfies 

U a (S u S 2 ,S 3 ) = U a {-8!, -8 2) -8 3 ), (23) 



and similarly for a — ► AB, one may verify directly that {{/^(fe), U B {k), V\(k), V 2 (k)} in Eq. (20i are all even under inversion 



of momentum, i.e. k — > (— fe). Physically, Eq. (23 1 expresses that the scattering potential for a set of lattice sites and the sites 
obtained upon a mirror transformation, as shown in Fig. [2] which is the case e.g. for a square lattice. In addition, one may verify 
that {[/^(fc), U B (k), Vi(fe), V^fc)} must all be real quantities for the same reason. 

Thus, we are finally able to give an analytical expression for the excitation energies for a two-component Bose-Einstein 
condensate with drag when taking into account arbitrary hopping and interaction between arbitrary sites. We find that 



1 r 



2[E 2 (k) + E\{k)\ + A[V 2 {k) - V 2 2 (k)} - 2[U 2 (k) + U 2 (k)} + 2a y/W) 



1/2 



(24) 



Ska - 2 

where we have introduced 

R{k) = 8[^ 2 (fe) + V 2 2 (k)][U A (k)U B (k) + E A (k)E B (k)]+4[V 2 2 (k) - V 2 {k)][U 2 A {k) + U 2 B (k) - E A (k) - E 2 B (k)} 

- 16Vi{k)V 2 {k)[E A {k)U B (k) + E B (k)U A (k)] + [E 2 A (k) + U 2 B (k) - E 2 B (k) - U 2 A {k)] 2 , (25) 



Eqs. (24 1 and ( |25| l represent one of our key results in this paper. Since there is no restriction on the sites involved in the hopping 
and interaction, the fc-dependence of the eigenvalues cannot be evaluated analytically in any straight-forward manner. However, 
the above closed analytical form for the excitation energies may serve as a basis for numerical investigations of the interaction 
between the two atomic species in the condensate. Below, we consider two limiting cases of particular relevance which allow 
further instructive analytical insight. 



B. Limiting case I: Nearest-neighbor hopping + on-site 
interactions 



We find that the terms in the Hamiltonian Eq. ( |T8| > may now 
be written as 



UqN 2 

2N„ 



+ N a (e 0a + T a ) 



UabNaNe 



(26) 



and we have introduced the basis vector 



t>k — [bkA, b-kA, b kB ,b_ kB ,bl A , b_ kA , b kB , 6l fcB ] T , 



(27) 
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FIG. 2: (Color online) The physical scenario expressed by Eq. t23 



where the 'T' superscript denotes the matrix transpose. The 
matrix A4k has an 8 x 8 structure, and reads 



Mk= 1 



Yk Xk 

upon defining the auxiliary matrices: 



(28) 



Yk = 








Fab 


° \ 





E£ 





Fab 


Fab 





Eg 





V 


Fab 





E B k) 


( 


F a 





F A b\ 


F A 





Fab 








Fab 





F B , 


\Fab 





F B 


/ 



(29) 



Upon introducing N a /N s — n a , we may write F AB = 
U AB\/ n A n B ■> F a = U a n a , and = t a J2s f 1 ~ e" lfe ' 6 ^) , 
Ek = e k + E a , a — A, B. By undertaking a diagonalization 
procedure, one obtains the excitation spectrum for the con- 
densed ground-state. Some care must be exercised in this pro- 
cedure, as the new quasiparticle operators in the diagonalized 
basis must also satisfy the boson commutation relations. As 
discussed previously, it is the matrix A^fccr 3 that must be diag- 



6 



onalized to obtain the quasiparticle excitation energies. Eval- 
uating the above determinant yields four distinct eigenvalues 
Afe = ±£ kcr , cr = ±1, where 



£ k.rr — 



e£(e£ + 2F A ) + eZ(e» + 2F B ) o 
2 2 



J[e£(e£ + 2F A )-e*(eB+2F B )]* 



1/2 



(30) 



Note that in the limit of two decoupled Bose-Einstein con- 
densates {Fab = 0) which are identical (Fa = F B = F, 
tA = ts = t), we regain the well-known single-component 
spectrum £ k = J e k (e k + 2F). The matrix D k now contains 
the excitation spectrum and reads (the choice of the order of 
the eigenvalues is arbitrary) 

D k = diag(d k ,d k ), d k = diag(£ fe+ , £ k _, -£ k+ , -£k-)- 

(31) 

Some comments are in order at this point. First of all, a 
similar approach to the condensed phase of a two-component 
Bose-Einstein condensate has been undertaken in both Ref. 
|[T3l and lfT2ll . However, the final answer for the diagonalized 
spectrum appears to be erroneous in Ref. lfT3l . where the ef- 
fect of the drag (interspecies coupling Uab) was completely 
disregarded in the excitation spectrum. Our results agree with 
the ones obtained in Ref. fl2l . The zero-temperature phase di- 
agram for a two-component Bose-Einstein condensate on an 
optical lattice was analytically constructed in Ref. [6]. More- 
over, it was pointed out in Ref. [ 14] that within the framework 
employed here (Bogolioubov approach) one is able to obtain 
the criteria that demarcates the transition from a superfiuid 
to Mott-insulator state, but one is not able to find the mani- 
festation of this phase transition in e.g. a sharp drop of the 
condensate fraction. 

We will now proceed to investigate the superfiuid velocity 
in more detail. The hydrodynamic flow in a Bose-Einstein 
condensate, and thus the superfiuid velocity, may be probed 
experimentally by stirring the condensate with, for instance, a 
blue-detuned laser beam as in Ref. 1151 . In the present case, 
we find two branches 1201 



= Vfcf . 



fc£fc<r|fc^0- 



(32) 



Below, we consider the one-dimensional case to obtain ana- 
lytically transparent results which should elucidate the basic 
physics. Straight-forward derivation leads to: 



aay(t A F A - t B F B y + AF 2 AB t A t B 

1 1/2 

■a 2 (t A F A + t B F B ) 



(33) 



This is consistent with the sound-like spectrum of Eq. ( 30 1 
in the long-wavelength limit k — ► 0. Note how the superfiuid 
velocity for each branch vanishes when the interaction param- 
eters U a , Uab m the problem are set to zero. Moreover, the 
superfiuid velocity v_ vanishes if one of the hopping matrix 



elements t A ort B vanishes, in which case the interspecies in- 
teraction parameter Uab is not relevant in the superfiuid ve- 
locity, such that v + reduces to the superfiuid velocity of a 
one-component Bose-Hubbard model. It is also interesting to 
generalize Eq. ( 30 1 to the case of particles moving in a con- 



tinuum, by substituting 



2m a 



(34) 



in which case the superfiuid velocity takes the form 




(35) 



Again, the result reduces to that of a one-component Bose- 
Hubbard model for the case where one of the species becomes 
immobile, i.e. either ttia or m B becomes infinite, and the su- 
perfiuid velocities vanish in the non-interacting case. In the 
continuum picture, we may also generalize Bogolioubovs ar- 
gument for the behavior of the excitations in the short- and 
long-wavelength limit. The limit of the long-wavelength lin- 
ear sound-like spectrum is roughly demarcated by a wavevec- 
tor which gives equal magnitude for the kinetic and potential 
energy terms in the quasiparticle dispersion relation. For com- 
ponent a, the crossover wavevector k c , a to the linear regime 
is given by 



2m n 



=>■ k C: a = \/2m a n a (U a + n & U a& ) ~ i. (36) 

where a denotes the other component in the condensate while 
£ a is the coherence length. The physical picture is then that 
the atoms of species a move as free particles on short length 
scales compared to £ a , while they move collectively at large 
length scales compared to £ a . Some other aspects of the su- 
perfiuid velocity for a two-component condensate with an en- 
ergy dispersion appropriate for the continuum were consid- 
ered in Ref. ifTBI . 

We now proceed to investigate in detail how the superfiuid 
velocity Eq. ( [33] l depends on the kinetic and potential energy 
terms in the problem. As seen, v a depends on the hopping 
parameters t a , the intraspecies interactions F a , and the inter- 
species interaction Fab- Upon choosing the parameters, we 
must ensure that the excitation energies remain real, as re- 
quired for a stable phase of two interacting atomic species. 



From Eq. ( 30 1, one infers that the solution may become imag- 
inary if the interaction j A b becomes sufficiently large. The 
criterion for a stable coexistent phase of the condensed phase 
for both atomic species reads O 



jAiB > Jab- 



(37) 



Let us first investigate how the two branches of the superfiuid 
velocity depend on the interspecies coupling. It is convenient 
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to rewrite Eq. p3\ in terms of dimensionless parameters as 
follows: 



TP 2 

t — + TP n 2 n— AB 

S,a — i'a r a<- L > P — tp tp 
r A?B 



Jab 

1A 1b' 



(38) 



It is interesting to note that for a fixed value of p, the tun- 
nel coupling amplitudes t a and the interaction parameters U a 
play the same role. The expression for the superfluid velocity 
remains the same under exchange of these two energy scales. 
The physical regime of the normalized interspecies coupling 
is now p € [0, 1], as demanded by Eq. (37 1. In Fig. [3] we show 
how the superfluid velocities in the two branches v± depend 
on the parameters in the problem. We give results for val- 



ues of p ranging from a weak interatomic scattering strength 
(p = 0.1) to a strong interaction (p = 0.9). As seen, the 
individual branches are not very sensitive to the value of p, 
but the two branches themselves differ qualitatively in their 
dependence on the hopping amplitudes and the potential en- 
ergy. In the case of two symmetric Bose-Einstein condensates 
(<^4 = (b) = £, one obtains from Eq. (38 1 that 



= [2<(1 + <Tyfp)\ 



1/2 



(39) 



The most interesting aspect of Fig. [3] is that the u_ branch 
vanishes as p — > 1. This means that a very small rotation or 
stirring of the condensate will trigger the a — — branch to 
become a normal fluid when p — > 1. 



f) = 0.1 



(9 = 0,6 



,3 = 0.9 




FIG. 3: (Color online) Plot of the superfluid velocity for the two quasiparticle branches v±. Above, we have used the dimensionless parameters 
("c = —t a F a a 2 and p = 7^ fl / (7a7s). The latter is a measure for the interaction between the two atomic components on the condensate. 



C. Limiting case II: Next nearest-neighbor hopping + off -site interactions 



We now go beyond the main approximation of Sec. |IIIB| and allow additionally for both next nearest-neighbor hopping and 
nearest-neighbor interactions. In this way, the interaction term in Eq. Q becomes: 



2 ^ S U t3kla( ]b\ a b) p b k pb la - [ Uafiblfilpbipbia + K 

ijkl a[3 a/3 i i^j 



j3 n ia n j(3 



(40) 



where U a p and U' a g are the on-site and nearest-neighbor in- 
teractions, respectively. In order to obtain transparent ana- 
lytical results, we consider the ID case, corresponding to a 
trapping potential which is elongated in a "cigar"-like shape. 



Proceeding in an equivalent manner as in the previous sec- 
tions, we finally obtain four distinct eigenvalues = ±£ka, 
a = ±1, which are identical to Eq. (BOb except that Fj — > F^, 
j € {A, B, AB}, where we have defined the kinetic energy 
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term 



2t a [l -cos(fca)] + 2t' a [l 



cos(2fcet)], a — A,B 

(41) 



and the potential energy terms 

= n a [U a + 2U' a cos(ka)], a ~ A,B 
F£ B = V^b~[Uab + 2U' AB cos(fco)]. 



(42) 



In Eq. (41 1, t a denotes the hopping parameter for nearest- 
neighbors while t' a denotes the hopping parameter for next 
nearest-neighbors. One now obtains the two branches of su- 
perfluid velocities which may be expressed through dimen- 
sionless quantities as follows: 



r, T = \l m A + m B + C7y/(^A - ^B) 2 + ^A^B, (43) 

where we have defined 

# a = (t a + K)(U a + 2U' a )n a a 2 , 



(u AB + 2u' AB y 

{U A + 2U' A ){U B + 2U' B y 



(44) 



Note that the above equations have exactly the same form as 
Eq. ([38]>, and that one obtains v — ► p, ^ a — > ( a in the limit 
{U' a , U AB } — > 0' as demanded by consistency. The stability 
condition for having a coexistent phase of the two superfluid 
branches is obtained by demanding that the eigenvalues are 
real, leading to the condition 



(4 



2i^)(e£ + 2F*)>A{f£ b ) 2 . 



(45) 



This is a generalization of the condition UaUb > U AB that 
arises from the standard assumption of only nearest-neighbor 
hopping and on-site interactions. Assuming e£ > 0, we may 
set = to find a more strict condition 



[U 



AB 



2U' AB cos(ka)}' 



■2U' A cos(ka)][U E 



2U' B cos(fca)] 



> 1 



(46) 



for the phase-coexistence regime. Thus, for a strong repulsive 
interaction between the atomic species A and B, one would 
expect that they do not coexist spatially but are instead sepa- 
rated into two distinct spatial regions. 



D. Microscopic parameters and drag between superfluid 
components 

We here derive explicit analytical expressions for the hop- 
ping and interaction parameters t a and U a p in our model. We 
will consider nearest-neighbor hopping and an optical lattice 
with intersite distance a. Our results are derived for the three- 
dimensional case, but are written down in a form which may 
be easily generalized to one or two dimensions. 

Starting from the definitions in Eq. Q, we obtain 




(47) 



which is consistent with Eq. (9) in Ref. (T2J. The definition 
of ujj a was given in Eq. pj. Above, j € {x,y,z}. Now, 
we present a derivation of the hopping term upon taking fully 
into account the trapping potential, which has been neglected 
in the literature so far. This is appropriate in a situation where 
the trapping potential has been turned off, allowing the con- 
densate to expand very slowly. Inserting the potentials and the 
Wannier functions into Eq. Q, we obtain 



*°=-n 

3 



dre -p a -R(o) 



2m a 

where we have defined 



V a • R(0) 



-p a -R(a) 



(48) 



3 3 

V a = (m a /2) + 2V 3 k 2 /m a )j. (49) 



After a shift of variables, we arrive at 



- (m Q 



,a 2 /4) 



j' 



(50) 



The ratio of the interaction term and the hopping term may 
now be evaluated straight-forwardly for any choice of micro- 
scopic parameters. Experimentally, it is possible to tune lat- 
tice parameters Vb and A through the laserlight setting up the 
optical potential. Defining the atom recoil energy 



E R , a = k 2 /2m a = 2n 2 /(m a \ 2 ) 

one may then define the tunable parameter 

Va 
Er,a 



(51) 



(52) 



which captures the effect of both the lattice well depth Vo and 
the lattice constant a = A/2. We simply denote Er.a = Er 
from now on. For later use, we note that for a cubic lattice and 
in the absence of a trapping potential, we obtain the relations 



Er 


= 3e" 37r V* 


U a 


a a 2 rtiA 


Er 


A 7r m a 


Uab 
Er 


_ a-AB 1 ,j 

A 7T 



(2TT^sm a /m A ) 3/2 , 
— -(1 + m A /m B )[ J 



3/2 



(53) 



upon choosing a positive sign for the hopping parameters. 
This fully determines the Bose-Hubbard parameters t a and 
U a for a given set of microscopic parameters. To access the 
physically allowed regime of qab, we define 



r, = ^^, V e[0,l). 

1A1B 



(54) 
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Symmetric case m Jm =1.0 




U A = U B 



Asymmetric case mym B =0.5 

0.3 r 




This is accomplished by virtue of our analytical expressions 
for both the quasiparticle energies (Sec. |III B\ and the param- 
eters in the effective Bose-Hubbard Hamiltonian derived pre- 
viously in this section. We here focus on the zero-temperature 
case, i.e. far away from the critical temperature, where our 
mean-field approach should be viable. 

We now derive an analytical expression for pd from the mi- 
croscopic Hamiltonian determined by Eqs. (18i, (26 1, (27i, 
(28j), (p9j. Our strategy is to let k 



k — m a v a in the Hamil 
tonian, leading to the Doppler-shifted energies 



e£ - m a v a ■ V fe £fc. 



(56) 



The energy eigenvalues may then be solved by expanding the 
characteristic polynomial in orders of v a , along the lines of 
IHn . At zero temperature, one obtains the following expres- 
sion for the drag coefficient: 



FIG. 4: (Color online) Plot of the hopping and interaction param- 
eters in the Bose-Hubbard model as a function of the trap depth 

s = V /E R . 



We plot in Fig. |4] the Bose-Hubbard parameters as a func- 
tion of the lattice well depth s for the case of symmetric 
{fix a/ Tn B = 1.0) and asymmetric (m^/raj = 0.5) two- 
component condensate. Moreover, a a /X = 10~ 3 , corre- 
sponding to a scattering length of a few nm for a typical ex- 
periment. As seen, the hopping amplitude becomes compa- 
rable to the interaction term only for optical lattice potentials 
Vq ~ En yielding s ~ 1 . We do not consider here very weak 
lattice potentials satisfying s <C 1, since the tight-binding 
model employed in the present paper no longer remains valid. 

As a practical application of our results for the excitation 
energies as well as the relation between the Bose-Hubbard pa- 
rameters and microscopic parameters, we now study the mag- 
nitude of the intercomponent drag coefficient pd in a uniform 
two-component BEC. In particular, we investigate what values 
Pd may realistically take for a relevant choice of microscopic 
parameters. The drag stems from a transfer of motion between 
the supercurrents for each component as a result of the inter- 
action jab, and vanishes in the case of two decoupled BECs. 
The free energy for a uniform two-component Bose-Einstein 
condensate may be written as [1 1 1 



_ AmAm B tAtB \ - 
Pd ~ N x N y N z a 



n B e£e»sm'(k x a) ^ 

£k,+£k,-(£k,+ + £k,-) ,: 



(57) 



F = F + V[p A v 



Pbv 2 b 



p d (v A -v B ) 2 }/2, (55) 



where Fq contains terms independent of the superfluid veloc- 
ities Vi for the two components and V is the volume of the 
system. The terms pj, j £ A,B represent the mass densities 
of the two components. 

In Ref. Q, an explicit expression was derived for the in- 
tercomponent drag for the case of small superfluid veloc- 
ities (much smaller than the critical ones) in the continuum 
limit, i.e. with free-boson dispersion relations e£ = k 2 /2m a . 
However, the drag between components on an optical lattice 
remains to be investigated. In what follows, we shall calculate 
Pd as a function of the microscopic parameters in the problem. 



Just like in the continuum limit treated in Ref. Q, we find 
that the drag coefficient is independent of the sign of the in- 
tercomponent scattering Fab- It is also seen that Eq. (57 i is 
always positive, pd > 0. Our results Eq. p7j ) may thus be 
considered as a generalization of the drag coefficient in Ref. 
Q to an optical lattice scenario. 

We now proceed to investigate the behavior of the drag co- 
efficient numerically on a 50 x 50 x 50 cubic lattice (Nj = 
50, Vj = Vq), j = {x, y, z}, which corresponds to the exper- 
imental setup of Ref. [ 19 1. Moreover, we fix n A — n B = V%, 
corresponding to an incommensurate filling as demanded for 
the superfluid phase. Let us define the normalized and di- 
mensionless drag coefficient pd/po, where po = tuaNa/V. 
For a fixed intracomponent interaction strength a a , the drag 
coefficient will thus depend on the strength of the intercom- 
ponent scattering rj, the mass ratio m^/raj, and the lattice 
well depth s. These microscopic parameters also determine 



the Bose-Hubbard parameters t a and U a /3 through Eq. (53 1. 
We present the dependence of Pd/po on mj/mg and s in Fig. 
|5]forboth a weak (rj — 0.2) and strong (r; = 0.8) intercompo- 
nent scattering. 

Qualitatively, it is seen that the plots are similar. With in- 
creasing lattice well depth s, the drag coefficient quickly di- 
minishes in size. It is interesting to note that the drag coeffi- 
cient is at its largest for a mass ratio m^/ms ~ 1, regardless 
of the value of s. This suggests that the velocity-drag effect 
between the components becomes most efficient when they 
have similar masses, which is reasonable. If pd ^ 0, it is pos- 
sible that the superfluid motion of one component induces a 
supercurrent in the other component purely by a drag effect, 
since the expressions for the supercurrents j a may be written 
as OH: 



3A = (PA - Pd)vA + PdVB, 3 B = (PB - Pd)v B + PdVA- 

(58) 

Thus, one may have j a ^ even with v a = 0. 
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drag has extremly important ramifications for the physics of 
these systems. Physically, a positive drag coefficient means, 



77 = 0.2 
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FIG. 5: (Color online) Contour plot of the normalized drag coeffi- 
cient Pd/po as a function of the mass ratio mA/rriB and the lattice 
well depth s. Two different values of r\ have been used, correspond- 
ing to weak (77 = 0.2) and strong (77 = 0.8) intercomponent scatter- 
ing. 



We end this section by briefly commenting on the positivity 
of pd that we find. In previous works on two-component Bose 
condensates, a negative drag has been found numerically in 
Monte Carlo computations |8 |, and starting from such a neg- 
ative value, highly unusual vortex states in a rotating multi- 
component Bose condensate, have been predicted which have 
no counterpart in the case pd > ifTOII . In particular, a su- 
perfluid bosonic density wave, corresponding to a vortex sys- 
tem which has the characteristics of a liquid and a solid at 
the same time, has been reported ifTUl . Therefore, a negative 



by virtue of Eq. 58 that a superfluid flow in one component 
of the condensate induces a co-directed, not counter-directed, 
flow in the other component. It should be noted that the neg- 
ative drag reported in Ref. [8] was obtained in a limit where 
the bosons on the optical lattice were very strongly interact- 
ing (essentially the hard-core boson limit) and the system was 
close to half-filling. Under such circumstances, one may ex- 
pect a backflow of one species of bosons when a boson of one 
component hops from one lattice site to another site occupied 
by the other component. Our approach, on the other hand, 
is essentially a weak-coupling approach which cannot capture 
such physics, and this issue warrants further investigation us- 
ing more suitable strong-coupling approaches. 



IV. SUMMARY 

In conclusion, we have investigated the excitation 
spectrum, superfluid velocity, and inter-component drag- 
coefficient for a two-component Bose-Einstein condensate on 
an optical lattice. We have derived analytical expressions 
for the excitation energies for arbitrary hopping and interac- 
tion between sites in Eqs. (24 1 and (25 1. We have investi- 



gated the excitation spectrum, superfluid velocity, and phase- 
separation condition in more detail for two important limiting 
cases of the general expressions. The critical superfluid ve- 
locity may be probed experimentally by using e.g. a laser as a 
macroscopic object to stir the hydrodynamic flow in the two- 
component Bose-Einstein condensate. Moreover, we have de- 
rived an analytical expression for the drag coefficient between 
the components when the condensate resides on an optical lat- 
tice in a weak-coupling approximation and found that it is al- 
ways positive. This means that drag induces co-directed flows 
of superfluid components, and not counter-flows, in the weak- 
coupling limit. We find that the transfer of motion from one 
supercurrent to another becomes most efficient for a mass ra- 
tio close to one for the two components, regardless of the lat- 
tice well depth or the intercomponent scattering strength. 
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